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We show that the jet structure of gamma-ray bursts(GRBs) can be inves- 
tigated with the tail emission of the prompt GRB. The tail emission that we 
p j| consider is identified as a steep decay component of the early X-ray afterglow 

q observed by the X-Ray Telescope on board Swift. Using a Monte Carlo method, 

£j ! we derive for the first time the distribution of the decay index of the GRB tail 

emission for various jet models. The new definitions of the zero of time and the 
time interval of a fitting region are proposed. These definitions for fitting the 
light curve lead us to a unique definition of the decay index, which is useful to 
l/") \ investigate the structure of the GRB jet. We find that if the GRB jet has a 

core-envelope structure, the predicted distribution of the decay index of the tail 
has a wide scatter and multiple peaks, which cannot be seen for the case of the 
uniform- and the Gaussian jet. Therefore, the decay index distribution gives us 
t — ! information about the jet structure. Especially if we observe events whose decay 



O 



index is less than about 2, both the uniform- and the Gaussian jet models will 
be disfavored, according to our simulation study. 



X 

Subject headings: gamma-rays: burst — gamma-rays: theory 



1. Introduction 



Gamma-ray burst (GRB) jet structure, that is, the energy distributio n E(Q) in the ultra- 



relativistic collimated outflow, is at present not yet fully understood (IZhang fc Meszaros 
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2002bl ). There are many je t models proposed in addition to the s i mplest uniform-jet model: 
the po wer-law jet model (IRossi et al.l 120021 ; IZhang fc Meszarosl l2002al). the Gaussian jet 
model (IZhang et al.l 12004 ) , the annular jet model (lEichler fc Leyinsonl I2004T ) . the multiple 



emittin g subshell model 



Kumar &: Piranl |2000| ; iNakamural 120001 ) . the two-component jet 



model (IBerger et al.l 120031 ) . and so on. The jet structure may depend on the generation 



process of the jet and therefore may provide us important information abou t the central 



engine of t he GRB. For example, in the collapsar model for long GRBs (e.g., IZhang et al. 



20031 . 12004 ). the jet penetrates into and break s out of the progenitor star, resulting in the 
E{9) oc 6~ 2 profile (ILazzati fc Begelmanll2005l ). For the compact binary merger model for 
short GRBs, hydrodynamic simulations have shown that the resultin g jet tends to have a 
flat core surrounded by the power-law-like envelope (jAloy et al.ll2005l ). 



In the pie-Swift era, there were many attempts to constrain the GRB jet structure. 
Thanks to the HETE-2, st atistical properties of long GRBs, X-ray-rich GRBs, and X- 
ray flashes were obtained ([Sakamoto et al.l 120051 ). which were thought to constrain the 
jet models (ILamb et al.l 12004 ). These observat i onal r e sults constrain va r ious jet models , 



such as the uniform-jet model ( Yamazaki et al 



2004al : lLamb et all 1200,4 iDonaghvi 120 06). 



the m ultiple subshell model (IToma et al.l l2005bl). the Gauss i an jet model (IDai fc Zhang 
2005h . and so on. For BATSE long GRBs, lYonetoku et~aD (hood) derived the dist ribu- 



tion of the pseudo-opening angl e, inferred from the Ghirlanda (IGhirlanda et al.l 12004 ) and 



Yonetoku (lYonetoku et al.l 12004 ) relations, as f(9j)d9j oc (L 2 d9j, whic h is co mpatible with 



that predicted by the power-law jet model as discussed in iPerna et al.l (120031 ) (however, see 
Nakar et al.l 12004). Afterglo w properties are also expected to constrain the jet structure 



[e.g.. iGranot fc Kumarl 120031 ); however, energy redistribution effects prevent us from reach- 
ing a definite conclusion. Polariz ation measurements of optical afterglows bring us important 
information (ILazzati et al.l 12004 ). 



In the Swift era, rapid follow-up o bservation reveals prompt GRBs followed by a steep de- 
cay ph ase in the X-ray early afterglow ( Tagliaferri et al]|2005 ; Nousek et all200~6 ; O'Brien et al 
2006al ). In the most popular interpretations, the steep decay component is the tail emis- 
sion of the prompt GRB (the so called high- latitude emission), i. e., the internal s hock 



origin (IZhang et al.l 120061 ; lYamazaki et al.l 120061; iLiang et al 



though there are some other possibilities (e.g . . iKobavashi et al 



2006 



Dvks et al.l 12005). al- 



2007; Panaitescu et al 



2006 



Pe'er et all I2OO6I : ILazzati k BegelmarJ hood : bado et al.ll2006l ). Then, for the uniform-jet 
case, the predicted decay in dex is a = 1 — /3, where we use the convention F v oc T _C V 1+/3 
( jKumar fc Panaitescu! l2000l ). For power-law jet case (E(8) oc 6~ q ), the relation is modified 
to a = 1 — (3 + (q/2). However, these simple analytical relations cannot be directly compared 
with observations, because they are for the case in which the observer's line of sight is along 
the jet axis and because changing the zero of time, which potentially lies anywhere within 
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the epoch where we see the bright pulses, substantially alters the early decay slope. 



Recently, lYamazaki et al.l (120061 ) (Y06) investigated the tail emission of the prompt 
GRB, finding that the jet structure can be described and that the global decay slope is 
not so much affected by the local angular inhomogeneity as it is affected by the global 
energy distribution. They also argued that the structured jet model is preferable, because 
steepening GRB tail breaks appeared in some events. In this paper, we calculate for the first 
time the distribution of the decay index of the prompt tail emission for various jet models 
and find that the derived distributions can be distinguished from each other, so that the jet 
structure can be more directly constrained than previous arguments. This paper is organized 
as follows. We describe our model in § [2j In § [3j we investigate the distribution of the decay 
index of the prompt GRB emission. Section @] is devoted to discussions. 



2. Tail Part of the Prompt GRB Emission 



We consider the sam e model as discussed in the previous works (Y06; lYamazaki et al. 



2004bl : iToma et alJ[2005aB l . The whole GRB jet, whose opening half-angle is A^tot, consists 



of iVtot emitting subshells. We introduce the spherical coordinate system (r, (p, t) in the 
central engine frame, where the origin is at the central engine and •d = is the axis of the 
whole jet. Each emitting subshell departs at time £^ p (0 < t^L < t dur , where j — 1, • • • , N tot , 
and tdur is the active time of the central engine) from the central engine in the direction of 
n (j) = <f U) ), and emits high-energy photons, generating a single pulse as observed. The 
direction of the observer is denoted by n Q b s = C#obs, V^obs)- The observed flux from the jth 
subshell is calculated when the following parameters are determined: the viewing angle of the 
subshell 9^ = cos -1 (n b S -n (j) ), the angular radius of the emitting shell A0^ b , the departure 
time tfe p , the Lorentz factor 7 (i) = (1 — Z^) -1 / 2 , the emitting radius rj$\ the low- and high- 
energy photon in dices and j3g\ the break frequency in the shell comoving frame z/q (j) 



( IBand et al.lll993l ). the normalization constant of the emissivity A u \ and the source redshift 
z. The observer time T = is chosen as the time of arrival at the observer of a photon 
emitted at the origin r = at t = 0. Then, at the observer, the starting and ending times 
of the jth subshell emission are given by 

-CO . 

<i) » , +0) i / i i _,2 nU) 2 



^ - C + ^J;(i + 7;^ )2 ) , (2) 

where 9^ = e®+A9%> b , 9 { i ] = max{0, 9^-A9^ h }, and we use the formulas f3 U) ~ 1-1/2^ 
and cos^ ~ 1 — 9 2 /2 for 7 (j) ^> 1 and 9 <C 1, respectively. The whole light curve from the 
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GRB jet is produced by the superposition of the subshell emission. 

Y06 discussed some kinematical properties of prompt GRBs in our model and found 
that each emitting subshell with 6^ 3) 3> produces a single, smooth, long- duration, dim, 

and soft pulse, and that such pulses overlap with each other and make the tail emission of 
the prompt GRB. Local inhomogeneities in the model are almost averaged during the tail 
emission phase, and the decay index of the tail is determined by the global jet structure, 
that is the mean angular distribution of the emitting subshell because in this paper all 
subshells are assumed to have the same properties unless otherwise stated. Therefore, we 
are essentially studying the tail emission from the usual continuous jets at once, i.e., from 
uniform- or power-law jets with no local inhomogeneity. In the following, we study various 
energy distributions of the GRB jet through the change of the angular distribution of the 
emitting subshell. 



3. Decay Index of the Prompt Tail Emission 



In this section, we perform Monte Carlo simulations in order to investigate the jet 
structure by calculating the statistical properties of the decay index of the tail emission. 
For a fixed-jet model, we randomly generate 10 4 observers with their line of sights (LOSs) 
^obs = (^obs, V^obs)- For each LOS, the light curve, F(T) of the prompt GRB tail in the 
15-25 keV band is calculated, and the decay index is determined. The adopted observation 
band is the low-energy end of the Burst Alert Telescope(BAT) detector and near the high- 
energy end of the X-Ray Telescope (XRT) on Swift. Hence, one can observationally obtain 
continuous light curves, beginning w ith the prompt GRB phase to the subsequent early 
afterglow phase (ISakamoto et al.ll2007l ). so that it is convenient for us to compare theoretical 
results with observations. However, our actual calculations have shown that our conclusion is 
not qualitatively altered, even if the observation band is changed, for example, to 0.5-10 keV, 
as usually considered for other references. 

For each light curve, the decay index is calculated by fitting F(T) with a single power- 
law form, cx (T — T*)~ a , as in the following (see Fig. [I ]). The decay index a depends on the 



choice of T* ( Zhang et al.ll2006l ; lYamazaki et alJbOQq jH Let T s and T e be the start and end 
time, respectively, of the prompt GRB, i.e., 



T s = min{T s ( 4j 



(3) 



1 Recently, iKobavashi fc Zhana (|2007l ) have discussed the way to choose the time zero. According to 
their arguments, the time zero is near the rising epoch of the last bright pulse in the prompt GRB phase. 
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maX { T end} 



(4) 



Then, we take T* as the time until 99% of the total fluence, which is defined by Stotai — 
f^" F(T') dT', is radiated, that is, 



Then, the prompt GRB is in the main emission phase for T < T*, while it is in the tail emis- 
sion phase for T > T*. The time interval [T a ,T&], in which the decay index a is determined 
assuming the form F(T) oc (T — T*)~ a , is taken to satisfy 



where we adopt q a = 1 x 10 2 and q b = 1 x 10 3 , unless otherwise stated. We find that in 
this epoch the assumed fitting form gives a well approximation. 

At first, we consider the uniform-jet case, in which the number of subshells per unit 
solid angle is approximately given by dN/dfl = N tot /(irA9% ot ) for i? < A0 tot , where A9 tot = 
0.25 rad is adopted. The departure time of each subshell is assumed to be homogeneously 
random between t — and t = tdur = 20 sec. The central engine is assumed to produce 
iVtot — 1000 subshells. In this section, we assume that all subshells have the same values of 
the following fiducial parameters: A9 suh = 0.02 rad, 7 = 100, r = 6.0 x 10 14 cm, a B = —1.0, 
(3 B = —2.3, hu' Q = 5 keV, and A = constant. Our assumption of constant A is justified 
as follows. Note that the case in which N subshells that have the same brightness A are 
launched into the same direction, but a different departure time, is equivalent to the case of 
one subshell emission with the brightness of NA. This is because in the tail emission phase, 
the second terms in the r.h.s. of Eqs. (1) and (2) dominate the first terms, so that the time 
difference effect, which arises from the difference of t^ep for each subshell, can be obscured. 
Hence, giving the angular distribution of the emission energy is equivalent to giving the 
angular distribution of the subshells with constant A. Also, Y06 showed that to obtain a 
smooth, monotonic tail emission as observed by Swift, the subshell properties, hv'^ and/or 
A U) , cannot have wide scatter in the GRB jet. Therefore, we can expect, at least as the 
zeroth-order approximation, that the subshells have the same properties. 

The left panel of Fig. [2] shows the decay index a as a function of $ Q bs- For -# b s < A9 tot 
(on-axis case), a clusters around ~ 3. On the other hand, when $ b s > A9 tot (off-axis case), 
a rapidly increases with i? b s - The reason is as follows. If all subshells are seen sideways 
(that is, 9^ ^> A0g^ b for all j), the bright pulses in the main emission phase followed by the 
tail emission disappear because of the relativistic beaming effect, resulting in a smaller flux 
contrast between the main emission phase and the tail emission phase compared with the on- 
axis case. Then T* becomes larger. Furthermore, in the off-axis case, the tail emission decays 




(5) 



F(T atb ) = q a<b F(T t ) 



(6) 
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more slowly (\dF/dT\ is smaller) than in the on-axis case. Then both T a — T* an d Th — T* are 



larger for the off-axis case than for the on-axis case. As can be seen in Fig. 3 of IZhang et al. 



( 120061 ). the emission seems to decay rapidly, so that the decay index a becomes large. The 
left panel of Fig. [3] shows a as a function of the total fluence S to tai which is the sum of the 
fluxes in the time interval, [T s , T e }. In Fig. [3], both a and Stotai are determined observationally, 
so that our theoretical calculation can be directly compared with the observation. 

A more realistic model is the Gaussian jet model, in which the number of subshells per 
unit solid angle is approximately given by dN/dQ = C exp(— $ 2 /2$j?) for ^ d ^ A0 to t> 
where C = iVt t/27r# 2 [l — exp(— A6^ ot /2^)} is the normalization constant. We find only a 
slight difference between the results for the uniform- and the Gaussian jet models. Therefore, 
we do not show the results for the Gaussian jet case in this paper. 

Next, we consider the power-law distribution. In this case, the number of subshells per 
unit solid angle is approximately given by dN/dVL — C[l + (^/^ c ) 2 ] -1 for ^ $ ^ A0 to t, 
i.e., dN/dVL ra C for ^ d < $ c and dN/dtt « C{d/d c )~ 2 for tf c < d ^ A6 tot , where C = 
(-Wtot/fl"#c)[ m (l + (A^tot/^c) 2 )] -1 is the normalization constant and we adopt $ c = 0.02 rad 
and A#tot = 0.25 rad. The other parameters are the same as for the uniform-jet case. 

As can be seen in the right panels of Figs. [2] and [3], both the ^obs^a and S to tai~tt diagrams 
are complicated compared with the uniform-jet case. When $ ohs < $ c , the observer's LOS is 
near the whole jet axis. Compared with the uniform-jet case, a is larger, because the power- 
law jet is dimmer in the outer region, i.e., emitting subshells are sparsely distributed near 
the periphery of the whole jet (see also the solid lines of Figs. 1 and 3 of Y06). If $ ODS 3> $ c , 
the scatter of a is large. Some bursts have an especially small a of around 2. This comes 
from the fact that the power-law jet has a core region (0 < $ < $ c ), where emitting subshells 
densely distributed compared with the outer region. The core generates the light-curve break 
in the tail emission phase, as can be seen in Fig. H] (Y06). In the epoch before the photons 
emitted by the core arrive at the observer (e.g., T — T s < 7.5 x 10 2 s for the solid line in 
Fig. @J, the number of subshells that contribute to the flux at time T, iV SUD (T), increases 
with T more rapidly than for the uniform-jet case. Then, the light curve shows a gradual 
decay. If the fitting region [T a , Tb] lies in this epoch, the decay index a is around 2. In the 
epoch after the photons arising from the core are observed (e.g., T — T s > 7.5 x 10 2 s for the 
solid line in Fig. H]), the subshell emission with 9 { v j) > $ obs + $ c is observed. Then N suh (T) 
rapidly decreases with T, and the observed flux suddenly drops. If the interval [T a , T b ] lies 
in this epoch, the decay index becomes larger than 4. 

To compare the two cases considered above more clearly, we derive the distribution of 
the decay index a. Here we consider the events whose peak fluxes are larger than 10~ 4 times 
of the largest one in all simulated events, because the events with small peak fluxes are not 
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observed. Fig. [5] shows the result. For the uniform-jet case (solid line), a clusters around 3, 
while for the power-law jet case (dotted line), the distribution is broad (1 < a < 7) and has 
multiple peaks. 

So far, we have considered the fiducial parameters. In the following, we discuss the 
dependence on parameters, r , 7, (3b, idur, and A# to t (It is found that the a-distribution 
hardly depends on the value of a^, A# su b, and u' within reasonable parameter ranges). At 
first, we consider the case in which r = 1.0 x 10 14 cm is adopted, with other parameters 
being fiducial. Fig. O shows the result. The shape of the a-distribution is almost the 
same as that for the fiducial parameters, in both the uniform- and the power-law jet cases. 
This comes from the fact that in a tail emission phase, the light curve for a given ro is 
approximately written as F(T; r$) « g(cT/r ), where a function g determines the light-curve 
shape of the tail emission for other given parameters. Then, the light curves in the case of 
r o = r o,i and r = r ,2, namely, F(T;r ,i) and F(T;r 0t2 ), satisfy the relation F(T;r 0t2 ) ~ 
F((ro^i/rQ^)T; ro,i). This can be seen, for example, by comparing the solid line with the 
dotted one in Fig. HI Hence, T*, T a , and Tb are approximately proportional to r ; in this 
simple scaling, one can easily find that a remains unchanged for different values of ro. 

Second, we consider the case of 7 = 200 and ro = 2.4 x 10 15 cm, with other parameters 
being fiducial. In this case, the angular spreading timescale (oc ro/7 2 ) is the same as in the 
fiducial case, so that the tail emissions still show smooth light curves, although the whole 
emission ends later, according to the scaling T e cx r (see the dot-dashed line in Fig. @J. 
Fig. [7| shows the result. For large 7, the relativistic beaming effect is more significant, so 
that the events in f^bs > A0 tot , which cause large a, are dim compared with the small-7 
case. Such events cannot be observed. For the power-law jet case, therefore, the number of 
large-a events becomes small, although the distribution is still broad (1 < a < 4) and has 
two peaks. On the other hand, for the uniform-jet case, the distribution of the decay index 
a is almost the same as for the fiducial parameter set, because the value of the decay index 
a in the case of ^ b s < A# to t is almost the same as that in the case of ^obs > A# to t- 

Third, we change the value of the high-energy photon index (3b from —2.3 to —5.0, with 
other parameters being fiducial. Fig. [H] shows the result. For the uniform-jet case, the mean 
value is (a) ~ 4, while (a) ~ 3 for the fiducial parameters, so that the decay index define d 



in this paper does not obey the well-known formula a = 1 — (3b ( jKumar fc Panaitescull2000l ). 
For the power-law jet case, the whole distribution shifts toward the higher value, and the 
ratio of the two peaks changes. In the tail emission phase, the spectral peak energy -E pea k is 
below 15 keV (see also Y06), so that the steeper the spectral slope of the high-energy side 
of the Band function, the more rapidly the emission decays, resulting in the dimmer tail 
emission (see the dashed line in Fig. H]). Then, the fitting region [T a , T&] shifts toward earlier 
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epochs, because T* becomes small. Therefore, the number of events with small a increases, 
and the number of events with large a decreases. Furthermore, we comment on the case in 
which fig is varied for each event in order to more directly compare with the observation. 
Here we randomly distribute /3b according to the Gaussian distribution with a mean of —2.3 
and a variance of 0.4. It is found that the results are not qualitatively changed. 

Next, we change the value of the duration time t^nr from 20 sec to 200 sec, with other 
parameters being fiducial. The epoch of the bright pulses in the main emission phase becomes 
longer than that in i<j ur — 20 sec. However, the behavior of the tail emission does not depend 
on tdur very much (see Fig. [9]) . Therefore, the distribution of the decay index a is almost 
the same as that for the fiducial parameters for both the uniform-jet case and the power- 
law jet case. Even if we consider the case in which t^ur is randomly distributed for each 
event according to the lognormal distribution with an average of log(20 s) and a logarithmic 
variance of 0.6, the results are not significantly changed. 

Finally, we discuss the dependence on A9 to t- Only the uniform-jet case is considered, 
because the structured jet is usually quasi-universal and because we focus our attention 
on the behavior of the uniform-jet model. The dotted line in Fig. [TD] shows the result for 
constant A6 to t =0.1 rad with other parameters being fiducial. We can see many events with 
large a. The large a is observed because for small A9 to t, although the off-axis events (i.e., 
A9 tot < -^obs) are still dim because of the relativistic beaming effect, a fraction of such events 
survives the flux threshold condition and are observable. Such events have large a > 5 (see 
the 4th paragraph of this section, which explains the left panel of Fig. [2]). This does not 
occur in the large-A# to t case. However, we still find in this case that there are no events 
with a < 2. We consider another case in which A#tot is variable. Here we generate events 
whose A0tot distributes as f A6tot d(A6 tot ) oc A6 tot - 2 d(A6 tot ) (0.05 < A6 tot < 0.4). Then for 
a given A9 tot , the quantities u' and A are determined by hu' = (A6> tot /0.13) -3 ' 6 keV and 
A oc (A^tot) -73 , respectively. Other parameters are fid ucial. If the model parameters are 



chose n in this way, the Amati and Ghirlanda relations (lAmati et al.ll2002l ; iGhirlanda et al. 



2004 ) are satisfi ed, and the eve nt rates of long GRBs, X-ray-rich GRBs and X-ray flashes 



become similar (iDonaghyl 120061 ) . The solid line in Fig. [TU1 shows the result. Again we find 



that there are no events with a < 2. 

In summary, when we adopt model parameters within reasonable ranges, the decay 
index becomes larger than ~ 2 for the uniform- and the Gaussian jet cases, while a significant 
fraction of events with a < 2 is expected for the power-law jet case. Therefore, if a non- 
negligible number of events with a < 2 are observed, both the uniform- and the Gaussian 
jet models will be disfavored. Furthermore, if we observationally derive the a-distribution, 
the structure of GRB jets will be more precisely determined. 
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4. Discussion 

We have calculated the distribution of the decay index, a, for the uniform-, Gaussian, 
and the power-law jet cases. For the uniform-jet case, a becomes larger than ~ 2, and its 
distribution has a single peak. The Gaussian jet model predicts almost the same results as 
the uniform-jet model. On the other hand, for the power-law jet case, a ranges between 
~ 1 and ~ 7, and its distribution has multiple peaks. Therefore, we can determine the jet 
structure of GRBs by analyzing a lot of early X-ray data showing a steep decay component 
that is identified as a prompt GRB tail emission. However, one of the big challenges in the 
Swift data for calculating the decay index in our definition is to derive the composite light 
curve of BAT and XRT. Since the observed energy bands of BAT and XRT do not overlap, 
we are forced to extrapolate one of the data sets to plot the light curve in a given energy 
band. To derive the composite light curve unambiguously for a prompt and an early X-ray 
emission, we need an observation of a prompt emission by current instruments, which overlap 
the energy range of XRT. 

The tail behavior with a < 2 does not appear in the uniform- and the Gaussian jet 
models; hence, it is important to constrain the jet structure. However, in practical observa- 
tions, such gradually decaying prompt tail emission might be misidentified with the external 
shock component, as expected in the pre-Swift era. Actually, some events have shown such 
a gradual decay, without the steep and the shallow decay phases, a nd their temporal an d 



spectral indices are consistent with a classical afterglow interpretation (lO'Brien et al.ll2006bl ). 
Hence, in order to distinguish the prompt tail emission from the external shock component 
at a time interval [T a , 7&], one should study the spectral evolution and/or the continuity 
and smoothness of the light curve (i.e., whether breaks appear or not) over the entire burst 
emission. 

In this paper, we adopt q a = 1 x 1CT 2 and g& = 1 X when the fitting epoch [T a , Tb] is 
determined [see Eq. (jBJ)]. Then, the prompt tail emission in this time interval is so dim that 
it may often be obscured by the external shock component, causing a subsequent shallow 
decay phase of the X-ray afterglow. One possible way to resolve this problem is to adopt 
larger values of q a and <#,, e.g., q a = 1/30 and qb = 1 x 1CT 2 , in which the interval [T a ,T&] 
shifts toward earlier epochs, so that the flux then is almost always dominated by the prompt 
tail emission. We have calculated the decay index distribution for this case (q a = 1/30 
and qb = 1 x 10~ 2 ) and have found that the differences between uniform- and power-law 
jets still arises as can be seen in the case of q a = 1 x 10~ 2 and qb = 1 x 10~ 3 , so that our 
conclusion remains unchanged. However, the duration of the interval, Tb — T a , becomes short, 
which might prevent us from observationally fixing the decay index at high significance. If 
Qa ^ 1/30, the emission at [T a ,T b ] is dominated by the last brightest pulse. Then, the light- 
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curve shape at [T a , T&] does not reflect the global jet structure, but reflects the properties of 
the emitting subshell causing the last brightest pulse. Another way to resolve the problem 
is to remove the shallow decay component. For this purpose, the origin of the shallow decay 
phase should be clarified in order to extract the dim prompt tail emission exactly. The other 
problem is contamination of X-ray flares, whose contribution has to be removed in order to 
investigate the tail emission component. In any case, if the GRB occurs in an extremely 
low-density region (a so-called naked GRB), where the external shock emission is expected 
to be undetectable, our method may be a powerful tool to investigate the GRB jet structure. 
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Fig. 1. — Example of how the decay index a is determined by the calculated light curve 
F(T). The start and end time of the burst are denoted by T s and T e , respectively. The time 
T* is determined by Eq. ([5]). The decay index a is determined by fitting F(T) oc (T — T*)~ a 
in the time interval [T a , T&]. 
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Fig. 2. — Decay index a as a function of f? bs; the angle between the whole jet axis and the 
observers' lines of sight. Red and green points represent events whose peak fluxes are larger 
and smaller than 10 -4 times the largest one in all simulated events, respectively. Left and 
right panels are for the uniform- and power-law jet cases, respectively. 
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Fig. 3. — Decay index a as a function of the total fluence Stotai, the sum °f the fluxes in 
the time interval [T s , T e ]. Red and green points represent events whose peak fluxes are larger 
and smaller than 10 -4 times the largest one in all simulated events, respectively. Left and 
right panels are for the uniform- and power-law jet cases, respectively. 
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Fig. 4.— Examples of light curves of the prompt tail emission in the 15-25 keV band 
for the power- law jet case and -# b s > $ c ("#obs = 0.27 rad and d c = 0.02 rad). The solid 
line shows the fiducial parameters. A bump caused by the core emission can be seen at 
T — T s ~ 7.5 x 10 2 s. The dotted, dot-dashed, and dashed lines are for r = 1.0 x 10 14 cm; 
tq = 2.4 x 10 15 cm and 7 = 200; and (3b = —5, respectively, with other parameters being 
fiducial. Time intervals [T a , T&] for each case are denoted by the thick solid lines. The flux 
is normalized by the peak value. 
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Fig. 5. — Distributions of the decay index a for uniform-jet (dN/dQ = const. ;solid line) and 
power-law jet (dN/dfl oc [1 + ($/ ^ c ) 2 ] -1 ; dotted line) models, respectively. We assume that 
all subshells have the same values of the following fiducial parameters: A0 su b = 0.02 rad, 
7 = 100, r = 6.0 x 10 14 cm, a B = -1.0, (3 B = -2.3, hu' = 5 keV, and A = constant. We 
consider events whose peak fluxes are larger than 10~ 4 times the largest one in all simulated 
events (red points in Fig. [2]). 
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Fig. 6. — Same as Fig. El but for tq = 1.0 x 10 14 cm. 
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Fig. 7. — Same as Fig. El but for 7 = 200 and tq = 2.4 x 10 15 cm. 
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Fig. 8. — Same as Fig. El but for (3b = 



-5.0. 
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Fig. 9. — Same as Fig. El but for tdur = 200 sec. 
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Fig. 10. — Distribution of the decay index a for the uniform-jet profile. The dotted line is for 
A# to t = 0.1 rad with other fiducial parameters. The solid line is for the variable- A# tot case, 
in which we generate events whose A^tot distributes as f& etot d(/\6 tot ) oc A# tot ~ 2 d(A# tot ) 
(0.05 < A9 tot < 0.4), and for given A0 tot , the quantities u' and A are determined by 
hv' = (A^tot/0.13)" 3 - 6 keV and A oc (A# to t) -7 ' 3 , respectively. Other parameters are fiducial. 



